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Abstract:A stand-alone efficient robust photovoltaic (PV)- 
Electrolyzer Scheme using a PV-Array, a DC-DC buck-boost 
converter, an Electrolyzer, a hydrogen tank, and an additional 
Fuel Cell (FC) stack feeding an standalone autonomous load 
is investigated and validated for efficient energy utilization, 
voltage stabilization and maximum power point tracking 
(MPPT) of the PV array under different insolation and 
temperature variations. The DC-DC Converter by controlled 
to ensure stable regulated operation. The PV generated power 
is used to electrolyze water to Hydrogen H 2 and Oxygen 0 2 , 
where hydrogen is stored in a suitably sized tank to be used 
later for the PEM-Proton Exchange membrane fuel cell 
(PEMFC). The PEMFC also acts as supplementary energy 
storage and load balancer to supply dynamic variable load 
whose value can vary from a minimum of 300W to a 
maximum of lKw. The Hydrogen H 2 flow to the FC depends 
on the dynamic load variations. The controller is fully 
regulated for efficient energy utilization as well as voltage 
regulation. The validated scheme components are modeled 
using Matlab/Simulink and fully integrated, to assess the 
unified system dynamic performance under varying 
conditions. 

Keywords — PV, Buck-Boost converter, Electrolyzer, and 
PEMFC. 

I. Introduction 

The ever-increasing demand for with depleted resources and 
massive energy consumption in developing countries created a 
dilemma of rising cost and limited supplies with soaring fossil 
fuel cost and global environment impact. This created and 
accelerated increased interest in renewable and alternative 
green enery/power generation systems. Photovoltaic (PV) 
power generation is the most promising renewable power 
generation technologies. The growth of PV power generation 
systems has exceeded the most optimistic estimation [1-2]. A 
significant problem for renewable energy generation system is 
need for energy storage for supplying the load demand during 
absence of Sun. Energy storage can be done via batteries such 
as Cr-Fe, Cd-Ni, lead-acid and vanadium redox flow batteries 
(VRB), [3-4]. The VRB is a promising green solution in 
practical applications [5], due to the advantages of large scale, 
simple structure, long life cycle and high SOC. After many 
technological advances, proton exchange membrane FC 
technology has now reached the test and demonstration phase. 
FC also shows great potential as a green power source, in the 
future because of many merits they have such as high 
efficiency, zero or low emission of pollutant gases, flexible 
modular structure and the rapid progress in FC technologies. 
However, each of the aforementioned technologies has its own 
drawbacks. For instance, wind and solar power are highly 
dependent on climate while FC needs hydrogen-rich fuel. 
Nevertheless, different alternative energy sources can 
complement each other to some extent. Multisource hybrid 



alternative energy systems (with proper control) have great 
potential to provide higher quality and more reliable power to 
customers than a system based on a single resource, [6] . 
The dynamic nonlinear Voltage-versus-Current (V-vs-I) 
characteristic of a PV array, there is a unique maximum power 
point (MPP) at which the PV operates at its maximum output 
power. However, the MPP moves under different atmospheric 
conditions because of the output voltage and current of a PV 
array vary as the radiation and temperature change. Operating 
the PV array at its maximum available power increases its 
efficiency resulting in a decrease in the operating costs. Since, 
PV arrays are still considered rather expensive compared with 
the utility fossil fuel generated electricity costs. The PV system 
owner desires to get maximum power from the PV array and to 
operate it at its highest conversion efficiency, [7]. In this 
paper control and management a system composed of a PV 
array, buck-boost converter, an electrolyzer, a H 2 tank, a 
proton exchange membrane fuel cell (PEMFC), and a lKw 
variable load is done. Each system component is modeled and 
integrated to demonstrate the system performance at different 
Insolation levels and load demand. A dynamic controller is 
designed and validated for Efficient Energy Utilization, 
Voltage Stabilization as well as load demand management. 

II. System Modeling 
1. PV Array Analysis 

The solar cell equivalent circuit is demonstrated by Fig. 1. The 
model (V-vs-I) characteristic equation describing the PV cell 
operation is given as, [8]: 




Fig. 1, PV- Cell Equivalent Circuit. 

Assuming: 1- The temperature is constant at 300° K (i.e. 27° C). 

2- The shunt resistance is too high, and 
consequently I sh is neglected. 

I D =I Q {exp[ A pv (V + IR Spv )]-l} 

(1) 

_ qV 



where 



PV 



AKT 
q: electron charge 
K: Boltzmann's constant 
A: completion factor 
Numerical methods have been used to calculate the cell 
parameters by using the experimental data at three different 
points on the I-V characteristics, namely, the short circuit 
current (0, I sc ), the open circuit voltage (V oc ,0) and the MPP 
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(Vp, Ip). Also dynamic resistances R S o~ "(— ) anc * ^shsc 

=-(^\ are obtained at the respective points. The equation 
V dl J I=hc 

of the generated current, which consists of N s cells in series 
and Np cells in parallel, is given by: 

Ir 



: Ip hg -Iog{exp[A g (Vp V +I pv R S g)]-l} 



R 



shg 



(2) 

The model of PV array in Matlab/Simulink software is shown 
in Fig. 2. The voltage-current characteristic curves of PV 
model under different irradiances are given in Fig. 3. 
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Fig. 2, PV array Matlab/Simulink model for varying Insolation 
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Fig. 3, V-I-P Relations between voltage, current and power of the PV 
array at different Insolation. 

The voltage at MPP is used as a reference value in the 
controllers in order to be tracked. The relation between the 
voltage at MPP and Insolation (§) can be expressed as a second 
order polynomial given by: 
V ref = V P 

= -0.0018(|) 2 + 0.2709tj) + 136.4725 

This relation is shown in Fig. 4 



Volt 




Insolation level 

Fig. 4, V-vs-cp Relation between voltage at MPP and Insolation level 

cp. 

2. DC-DC Converter-Buck-Boost 

An ideal Buck-Boost converter is used to interface the PV 
array source to the electrolyzer Load Cell. With the 
assumption of ideal circuit elements, two switched models are 
shown in Fig. 5. State variables for this Buck-Boost converter 
are chosen as the inductor current, I L =Xi, and the capacitor 
voltage, V c = X 2 . State-space-averaged equations in matrix 
form are, [11]: 



U, = V s ,P S 



>- 



"' fw 



4f- 



Fig. 5, Switched models for the ideal Buck-Boost converter. 
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3. Electrolyzer 

The Electrolyzer is a device that produces hydrogen and 
oxygen from water. Water-Eelectrolysis can be considered as a 
reverse process of a hydrogen PEM- fuel cell. The electrolyzer 
converts the DC electrical energy into chemical energy stored 
in hydrogen. Currently, there are three principal types of 
electrolyzer available in the market: alkaline, PEM, and solid 
oxide, [12]. 

The alkaline and PEM-Electrolyzers are well established 
devices with thousands of units in operation, while the solid- 
oxide Electrolyzer is as yet unproven [13]. Alkaline water 
electrolysis is the dominating technology today. In this section, 
the principle of alkaline water electrolysis is reviewed. An 
alkaline Electrolyzer uses potassium hydroxide (KOH) as 
electrolyte solution for transferring hydroxy 1 ions. Fig. 6 shows 
the schematic diagram of an alkaline Electrolyzer. 
At the cathode, two water molecules are reduced 
electrochemically (by the two electrons from the cathode) to 
one molecule of H2 and two hydroxyl ions (OH ). The reaction 
at the cathode can be expressed as: 
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DC Voltage 




Fig. 6, Schematic diagram of an alkaline type Electrolyzer. 

2H 2 0 + 2e" => H 2 1 +20H (5) 
Under the external electric field, the hydroxyl ions are forced to 
move towards the anode through the porous diaphragm. At the 
anode, the two molecule hydroxyl ions, losing two electrons to 
the anode, are discharged into 1/2 molecule of 0 2 and one 
molecule of water. The chemical reaction at the cathode is 

20H - ^>Vi0 2 1 + H 2 0 + 2e" (6) 
By combining the above two equations, the overall chemical 
reaction inside an electrolyzer is 

2H 2 0 => H 2 1 +Vi 0 2 1 (7) 
The current density of alkaline-Electrolyzers is normally less 
than 0.4 A/cm 2 . The energy conversion efficiencies can range 
from 60-90%. Without auxiliary purification equipment, 
purities of 99.8% for H 2 can be achieved. Alkaline electrolysis 
technology can be implemented at a variety of scales from less 
than 1 kW to large industrial electrolyzer plant over 100 MW as 
long as there is a proper DC electricity supply [14]. The model 
of an electrolyzer stack developed for this study is based on the 
empirical /-V equation reported in [15], described as 



I + k elec ln (- 



-1 + 1) 

(8) 



where Velec,cell is the cell terminal voltage (in volts), 
Vrev is the reversible cell voltage, 
r\ (in ohms square-meter) and r2 (in ohms square-meter 
per degree Celsius) are the parameters for the ohmic 
resistance inside the electrolyzer, 
ke lec (in volts), k T i (in square meters per ampere), k T2 
(square-meter degrees Celsius per ampere), and k T3 
(square-meter degree Celsius square per ampere) are 
the parameters for the overvoltage, 
A is the area of the cell electrode (in square-meters), 
/ is the electrolyzer current (in amperes), 
and Tis the cell temperature (in degrees Celsius). 
The V-I characteristics of the electrolyzer model used in this 
study at different cell temperatures are given in Fig. 7. At a 
given current, the higher the operating temperature, the lower is 
the terminal voltage needed. 

According to Faraday's law, hydrogen production rate for an 
Electrolyzer cell is directly proportional to the electrical current 
in the equivalent Electrolyzer circuit [16]. 



2F 



(9) 



where i e is the electrolyzer current, 

n c is the number of electrolyzer cells in series, 

and r\ F is the Faraday efficiency. 
Assuming that the working temperature of the electrolyzer is 40 
°Q Faraday efficiency is expressed by [8] 



r, F =96.5e (a09/i - 75 - 5/ ^ (10) 

Using equations (9) and (10), a simple electrolyzer model is 
developed using Simulink. 




Fig. 7, V-I characteristics of the H2-Electrolyzer Cell model 
under different temperatures. 

4. Hydrogen Storage Tank (HST) Modeling 

Physical hydrogen storage is one of the storage techniques that 
use tank to store compressed hydrogen. The model of HST is 
described according to equation (11), [17] and using the ratio 
of hydrogen flow to the tank, directly obtained the tank 
pressure. 

P b -P bi =z((n H2 -W H2 )RT b /M H V b ) (11) 

Where: R is universal gas constant [J (kmol °K) _1 ] 
P b is Pressure of tank [PJ 
P bi is Initial pressure of the storage tank [PJ 
W H2 is hydrogen flow rate delivered to FC, 
M H2 is hydrogen mass 
T b is Tank operating temperature [°K] 
V b is volume of the tank [m 3 ] 
z is compressibility factor as a function of tank 

pressure 

5. Proton Exchange Membrane Fuel Cell (PEMFC) 

FC converts oxygen and hydrogen chemical energy to 
electrical energy. The amount of output voltage in each cell is 
small. Therefore, cells are connected in series to gain a higher 
output voltage. Three basic are as which affect the output 
voltage are ohmic, activation and concentration polarization 
[18]. FC output voltage is mentioned in polarization curve 
calculated by: 

VFC □ □ Ee rnst □ V act D □ YhmicH □ 

(12) 

The open circuit voltage in FC which is denoted by E mmst can 
be defined as follows: 

E nernst = 1.229-0.85 * 10 3 (T - 298.15) + 4.3085* 10" 5 T[ln(P H2 ) + 0.5 ln(P 02 )] 

(13) 

where:T is refers to FC operation temperature, [°K]; 

P H2 & P02 are partial pressure of hydrogen & oxygen, 

[atm]. 

Ohmic polarization shows ohmic voltage drop due to resistive 
losses in the cell which is caused by resistance against 
electrons and ions movement. The ohmic voltage drop can be 
defined as follows: 

Vohmic =1fc(Rm+Rc) 

(14) 

where: / FC is refers to FC current; 

R c is resistance of FC electrodes against electrons 
movement which is constant; 
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R M is FC membrane resistance against ions passing, 
which equal p M (l/A) , [19]. 
Activation polarization shows activation voltage drops in FC 
anode and cathode electrodes. These drops are caused by the 
slow reactions on the surface of electrodes which are calculated 
by: 

V act = -Ei + ^ 2 T + ^Tln(C 02 ) + ^Tln(i FC ) 

(15) 

where C0 2 is the oxygen concentration rate on the catalyst 
surface. 

Concentration polarization shows voltage drop due to reduction 
in density of reaction materials which is called mass transport 
losses. To determine voltage drop model due to mass transport, 
maximum current density must be gained (/ max ). Voltage drop 
corresponding to mass transport can be calculated by [20] : 
V con =-Bln(l-(J/J ma x)) 

(16) 

Considering all FCs voltage drops and for n cells which are 
connected in series by forming a stack, output power of FC 
stack is defined as follows: 
Pfc = Npc Vpc Ifc 

(17) 

where N FC represents the number of cells in series. 

The FC current / FC can be determined as a function of hydrogen 

flow rate to FC and can be obtained by [21]: 

ifc = (2F/JV fc )(Wh2/Mh2) (18) 
where: W H2 is hydrogen molar flow rate delivered to FC, 

M H2 is hydrogen molar mass and F refers to Faraday's 
constant [C/kmol]. 

III. Power Management 

An overall control strategy for power management among PV 
array with buck boost converter, electrolyzer and Hydrogen 
tanks PEMFC and variable load is designed. Fig. 8 shows the 
block diagram of the overall control strategy for the proposed 
system. PV electricity generation unit controlled for voltage 
regulation controller. 
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Fig. 8, Block diagram of the overall PV-FC-Electrolyzer Energy 
Scheme 

IV. Digital Simulation Results and Discussion 

System input data includes load electricity demand and solar 
radiation during 24 h. A load varying between 300W and lkW 
along 24 hours is assumed. The demand of load electricity 
during 24 h is depicted in Fig. 9. The rates of hydrogen 
production by EL, consumed hydrogen by PEMFC and stored 



hydrogen in HST are shown in Fig. 10. Figures 9 and 10 
demonstrate the hydrogen consumed by the distributed load 
during the night (from 12:00 AM till 6:30 AM). As the sun 
rises it is clear that the amount of stored hydrogen that is 
consumed by the load decreases, while the hydrogen level in 
the tank is increasing. This increase in hydrogen level 
continues until 5:00PM, after which as the insolation level 
decreases, the H 2 level in the tank starts decreasing while 
supplying the load demand. Fig. 11 shows the simulation 
system. Fig. 12 shows the insolation level with change from 
low to high at t=0.1 sec and change from high to low at t=0.15 
sec. While, Fig. 13 shows the simulation results of voltage and 
current of photovoltaic at change in insolation level. It shows 
the controller performance as described in equation (3). 
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Fig. 9, 24 Hour-Load-Demand. 




Fig 10, The Hydrogen flow 
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Fig. 11, Simulation of the system. 
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Fig. 12, Change of insolation. 



IV. Conclusion 

An efficient regulation and control system for Power/demand 
load regulation of the flow of power generated by PV array to 
supply a variable load demand is designed. To ensure 
continuous supply of load, storage system is designed. The H2- 
storage system consists of an Electrolyer, H 2 tank with PV 
source and Fuel Cell storage and backup energy fuel cell. 
Maximum power utilization using the PV-MPPT based voltage 
regulator is fully validated under varying condition. Unified 
System components are modeled, and dynamic simulation 
results validated the robustness and damping behaviour under 
varying insolation and load conditions. He proposed scheme 
can be extended to PV Farm Utilization and 
Village/Island/V illage Electricity. 
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Fig. 13, the simulation results of Voltage and Current 
photovoltaic at change of insolation. 

Fig. 14 shows the output voltage converter at MPP voltage and 
voltage regulator. It compares the converter output with MPP 
control, and when applying constant voltage control. The 
output voltage of the PV array is regulated at 145 V. This is 
because the variation of the output voltage, when following 
MPP, reaches 141V at a 20% Insolation level, and 146.5 at 
80% Insolation level. Hence regulating the output voltage at 
the 146.5V outweighs its slight variation by 2.81% only 
between the maximum and minimum Insolation level. Figures 
15 and 16 show the converter supplying the Electrolyer, 
revealing the rate of hydrogen flow to be stored in the tank. 
These figures show the time required for filling the H 2 tank is 
about 4 hours. 
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(b) Voltage Regulator 
Fig. 14, Output voltage of Converter 







































































































— r 




















/ 





















0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 
Time (sec) 

Fig. 15, Output of hydrogen flow from Electrolyer 
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Fig. 16, the hydrogen storage tank 
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